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Abstract 

In order to confirm the picture of domain-wall excitations in the hidden an- 
tiferromagnetic order of the Haldane phase, the structure of the low- lying 
excitations in the S = 1 antiferromagnetic Heisenberg chain is studied by a 
quantum Monte Carlo method. It is confirmed that there exists a finite en- 
ergy gap between the first- and the second-excited states at k = tt as well as 
between the ground state and the first-excited state at k = tt. In the thermo- 
dynamic limit, the second-excited state at k = tt is separated from the ground 
state by the gap which is three times as large as the Haldane gap. From the 
size dependences of the low-lying-excitation energies, the interactions between 
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the elementary excitations in the excited states are concluded to be repulsive. 
PACS numbers: 75.10.Jm, 05.30.-d, 75.40.Mg 
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The structure of low-lying excitations of the S = 1 antiferromagnetic Heisenberg chain 
is one of the main interests in the study on the Haldane system. []I],|2|] The valence-bond- 
solid (VBS) model introduced by Affleck, Kennedy, Lieb, and Tasaki (AKLT), which 
gave a clear-cut physical picture of the Haldane massive state, stimulated several authors 
to study the nature of the excited states. In order to describe the elementary excitations, 
Knabe || considered a triplet bond constructed from two spin- 1/2 degrees of freedom in the 
VBS background, which is now called a crackion, while Arovas, Auerbach, and Haldane || 
discussed a domain wall in the hidden antiferromagnetic order. Later, Fath and Solyom || 
demonstrated that both of the defects have a solitonic nature and actually give the same 
dispersion relation that well reproduces the elementary excitation spectrum of the AKLT 



model. PU On the other hand, proposing a new Monte Carlo technique, Takahashi flO 
pioneeringly calculated the lower edge of the excitation spectrum as a function of momentum 
k for the Heisenberg chain. There he further suggested that the lowest excitation of low 
momentum may be a scattering state of two elementary excitations of momentum k ~ ir. 
Motivated by his calculation and based on an idea of the hidden domain wall, Gomez- 



Santos [IT] made a variational approach to the Heisenberg chain and showed that the lowest 
excitations are single-particle-like in the vicinity of k = 7r, whereas two-particle-like near 
k = 0. Thus the lower edge of the excitation spectrum has almost fully been investigated 
so far and it is probable that the elementary excitations are more or less identified with the 
moving hidden domain wall. 

Although nowadays the higher excitation with arbitrary momentum is generally believed 
to be a scattering state of the elementary domain walls, quantitative investigation of them 
has not yet been performed well. Developing quite different numerical treatments, White and 



Huse (T2j and the present authors [jnj qualitatively pointed out that the lowest excitations 
as a function of k constitute an isolated band at large fc's, while they coincide with the lower 
edge of the two-domain- wall continuum at small fc's (Fig. [1]). Their suggestion motivates us 



to calculate the second-lowest excitations. Making use of the Lanczos method, Takahashi |TJ[ 
calculated S(q, u) of the Heisenberg chain of short length as a series of 5-function peaks at 



each excited state, which suggests that the lowest excitations are actually separated from 
the upper continuum at large g's. For the AKLT model he found that the excitation energy 
of the second-lowest state at k = 7r is three times as large as one of the lowest state at 
k = 7r. However, the chain length L < 20 he treated was not long enough to confirm the 
same scenario in the thermodynamic limit at the Heisenberg point where the correlation 
length is much longer than that of the AKLT model. Thus the expected multi-domain-wall 
excitation energies in the thermodynamic limit, 

E 2 (0)-E G = 2A, £ 2 (vr)-£ G = 3A, (1) 

were not definitely obtained. Here, Ei(k) is the / th eigenvalue in the k- momentum space, 
Eq = Ei(0) is the ground-state energy, and A = Ei(n) — E G is the Haldane gap. 

In order to study the structure of the low-lying energy levels, we exploited an efficient 
quantum Monte Carlo method. In principle the structure factor S(q, uo) can be obtained from 
the spin correlation function along the Trotter axis, C(r), which is called the imaginary time 
correlation function. However, due to an intrinsic difficulty in the numerical treatment of 
transforming the C(r) to S(q,u), certain statistical methods such as the maximum-entropy 
technique, []T3-ff7| have to be introduced for practical calculations at finite temperatures. On 



the other hand, the lower boundary of the spectrum is simply obtained from a logarithmic 
plot (log plot) of C(r). For the Haldane systems this simple approach works well and 



a precise value of the lowest state at each momentum has been obtained. [|l8H2ll Let us 
outline the log-plot method. 

We consider a periodic chain of L spins described by the Hamiltonian 

L 

H = Sj ■ Sj + i ; Sl+i = Si , (2) 

3=1 

where Sj is the S — 1 spin operator. Let us denote the / th eigenvector and eigenvalue of 
7i in the /c-momentum space by \l; k) and Ei(k) [E\{k) < E2(k) < ■ • •], respectively. When 
the Hamiltonian has the translational symmetry, i.e., [H,T] = 0, 

H\l]k) = Ei{k)\l]k), T \l; k) = e ik \l; k) , (3) 



where T is the translation operator. The dynamic structure factor in a real frequency do- 
main, S(q,u), is the Laplace transformation of the imaginary-time spin correlation function 
of S z q = L~ l J2f=i S-e [qj : 

S(q,r) = (e^S z q e-^S z _ q ), (4) 

where (A) = Tr[e~^A]/Tr[e - ^ w ] denotes the canonical average at a given temperature 
f3~ x = k^T. When the system has some conserved quantities, the Hamiltonian is block- 
diagonalized. Since the total magnetization, M = J2i S z , is a conserved quantity in the 
present system, S(q, r) is independently defined in each subspace with a given M. Using 
the complete vector set \l; k) in each subspace, S(q,r) is represented as 

E,„e-^« (l;k\S z q \l';k + q) * e -^(*+</)-£«(fc)] 
Siq > T) = — • (5) 

Thus S(q, t) as a function of r generally exhibits a complicated multi-exponential decay. At 

a sufficiently low temperature, S(q, r) is given as 

S(q,r) = £ |(1; k \S z q \l; k + q )\ 2 e-rm^-E^)] ? (g) 
I 

where ko is the momentum at which the lowest state in the subspace is located. Now it is 
reasonable to approximate Ei(k + q) — Ei(k ) by the slope — dln[S(q, r)]/dr in the large-r 
region satisfying 

1/1 • U.\ CzL. I, i «\|2 

(7) 

for an arbitrary n. When the excitations constitute a isolated band and its spectral weight 
| (1; k \S q |1; k + q)\ 2 is large, the inequality (0) is well justified in a wide region of r. In fact 
these conditions are satisfied in the Haldane systems. That is why the method did work 
quite well for the lowest excitations of the present system especially in the large-g region. [TJ[ 



l(i;fcol^ 


n; k + q)\ 2 


l(i; 


Qz 


l;k + qW 



Taking a large enough Trotter number, we might, in principle, extract the higher-lying 
levels from the log plots of C(r) for proper regions of r. However, such an attempt has turned 
out to be unfeasible with the present numerical facility. Thus we use the log-plot method 



for the lowest level at each g, but in the subspaces with various values of M. Combining 
the calculations in these subspaces, we can construct the low-energy structure and obtain 
the higher-lying levels. Figure [I] illustrates the probable spectrum of the low-lying states of 
the system, which is based on the qualitative arguments [ |T2| , |I"3f and the calculations for the 
short chains of L < 20. 

In Fig. H we show log plots of S(q, r) calculated in the subspaces of M — 0, M = 1, and 
M = 2 as a function of r at various values of q, where the lattice constant was set equal 
to unity. It seems that in general lnfS^g, r)] shows a better linearity at large g's than at 
small g's. We have performed at least a million Monte Carlo steps to obtain S(q,r). In 
most cases the temperature (/5J) _1 and the Trotter number n have been set equal to 0.02 
and 200, respectively. We have checked that the thus-obtained data are reliable enough to 
represent the ground state properties. The numerical precision of the raw data amounts to 
two digits or more. 



As has been mentioned in our previous work, |jl3[| a single spin flip create a domain wall 
of the hidden order parameter. Thus the lowest-excitation energy is the gap between the 
lowest singlet state (S = 0) and the lowest triplet state (S = 1). Because the S = state 
contains a level of M = 0, the gap immediately above the ground state is obtained as the 
energy difference between the two lowest states in the M = subspace except for the case of 
q = 0. For q — 0, S* commutes with the Hamiltonian and therefore S(q, r) does not depend 
on t. The excitation energies from the ground state for the M = states are obtained as 

AE (q) = A (g) (g^0). (8) 

Here AE m (q) is the energy difference between the lowest-excited state with a momentum q 
in the M = m subspace and the ground state, and A m (g) is the present numerical finding, 
— dln[S(q, r)]/dr, in the subspace of M = m under the condition (^). Takahashi showed in 
his recent work |T^| that the spectral weight of the lowest excitation at each q is extremely 



large especially at q > 0.37T where the lowest excitation is separated from the upper con- 
tinuum. |TTlJl^| This fact causes the fine linearity of ln[,S(g, r)] at large g's but the less fine 



linearity at small g's, as was observed in Fig. 0(a). In Fact, in our previous attempt [18 to 
obtain the lower edge of the spectrum, the data at small g's were not so satisfying as ones 
at large g's from the point of the numerical precision. We note that the inelastic-neutron- 
scattering measurements |23|] actually revealed the clear isolated band for the 5* = 1 Haldane 
material Ni(C 2 H 8 N 2 ) 2 N0 2 (C10 4 ). 

In the subspace of M = 1, the present approach brings us the energy differences between 
the lowest triplet state (k = 7r) and other triplet states with an arbitrary momentum. The 
excitation energies from the ground state for the M = 1 states are obtained as 

AE 1 (q) = A (7r) + A 1 (7T-g). (9) 

In principle, AE (q) obtained from the calculation for M = should agree with A£' 1 (g) 
obtained from one for M = 1. However, a certain difficulty in the present treatment prevents 
us from reaching the definite coincidence between them. In the M = subspace, the lowest 
state at each q has a relatively large spectral weight, while in the M = 1 subspace, this 
seems not to be the case. Actually, for M = 1, the effect of the multi-exponential decay 
more clearly appears in the log plot of S(q, r) especially at small g's. But, even in the case 
of M = 1, the log plot of S(q, r) still gives an almost straight line in the vicinity of q = it, 
which results in a precise estimate of AE (0) = AEi(0) = A (ir) + Ai(7r). 

In the subspace of M = 2, the present calculation brings us the energy differences 
between the lowest quintuplet state (k = 0) || and other quintuplet states with an arbitrary 
momentum. The excitation energies from the ground state for the M = 2 states are obtained 
as 

AE 2 (q) = AE (0) + A 2 (g) = A (tt) + A l {n) + A 2 (g) . (10) 

Here nonlinearity of ln[S*(g, r)] persists in a relatively wide range of r even at large g's, which 
prevents us obtaining a full dispersion curve of the second-excited states. 

We plot in Fig. [5] the lowest and the second-lowest eigenvalues, AEi(q) and A£' 2 (g), as 
a function of g for the L = 64 chain, which have been obtained using Ao(vr) and Ai(7r) at 



L = 64. Here all the errors arise in estimating the slope of ln[>S(g, r)] rather than come from 
the raw Monte Carlo data. The lowest eigenvalues AEi(k) were determined, using both of 
S(q, t)'s with M = and M — 1, so as to minimize the numerical ambiguity. Although 
AE2(q) has not successfully been obtained in the region of q ^5 0.77T, we here clearly confirm 
the existence of the isolated band and observe the lower edge of the three-domain-wall 
continuum as well as one of the two-domain-wall continuum. The overall behavior is almost 
the same as that for L = 128 except in the vicinity of q = and q = ir, and therefore we 
can believe that we are observing the bulk behavior. 

In Fig. |] we show size dependences of the first-excitation energy at the zone boundary, 
AE (tt) = A£ , 1 (7r), the first-excitation energy at the zone center, AEi(0) = AE 2 (0), and 
the second-excitation energy, AE 2 {j{), which are supposed to be the bottom of the single- 
domain-wall band, one of the two-domain-wall continuum, and one of the three-domain-wall 



continuum, respectively. Here the symbols x represent the Takahashi's data |14j] for AE^ty) 
obtained through a different method, which are somewhat inconsistent with our finding 
beyond the numerical uncertainty. We here observe that the relation ([I]) comes to hold as 
L increases. Therefore, we conclude that the low-lying excitations of the present model are 
regarded as the domain-wall excitations in the hidden antiferromagnetic order. 

The chain-length dependences of the energy of the domain-wall-scattering state suggest 
that there exists a long-range repulsive interaction between the domain walls in the excited 



states, which is qualitatively consistent with a variational calculation. |24| We note that the 



interaction between the domain walls is contrastingly attractive in the ground state. [25 In 



the AKLT model, ]14j the relation (|l|) almost holds even at L = 16. It is well known that 



the spin-spin correlation length is l/ln3 ~ 0.91 at the AKLT point, |4j] while it is estimated 
to be 6.2 at the Heisenberg point. ]2£| Thus the spatial extension of the domain wall is 
expected to be much larger in the Heisenberg model than in the AKLT model, to which 
the present significant size dependence is attributed. Smearing the domain wall over three 
lattice sites, Scharf and Mikeska || obtained a variational bound for the Haldane gap of the 
AKLT model which coincides with the exact-diagonalization result within 1% error. The 

8 



spatial extension of the domain wall in the Heisenberg model may reach more than twenty 
lattice sites. 
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FIGURES 

FIG. 1. Illustration of the spectrum of the low-lying excited states as a function of momentum 
for the infinite chain. 

FIG. 2. Logarithmic plots of S(q, r) versus r at various values of q for the L = 64 chain: (a) 
M = 0, (b) M = 1, and (c) M = 2. 

FIG. 3. The lowest and the second lowest eigenvalues as a function of q for the L = 64 chain, 
where 0> ^> an d a denote the results obtained from S(q,r) calculated under M = 0, 1, and 2, 
respectively. 

FIG. 4. Size dependences of the excitation energy of the first-excited state at the zone boundary, 
AEq(tt), one of the first-excited state at the zone center, AE\(0), and one of the second-excited 
state at the zone boundary, AE^tt). The black symbols represent the L — > oo extrapolated 
values. The symbols x represent the Takahashi's calculations [14] for A£^(y) obtained with an 
exact-diagonalization technique. 
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